arXiv:1501.04472vl [cond-mat.soft] 19 Jan 2015 


1 HISTORY DEPENDENCE IN GRANULAR MATTER 
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The mechanical properties of a granular sample can depend on the way the packing was prepared. However, the variables which 
store this information are often unknown. Here we present an X-ray tomography study of tetrahedra packings prepared with three 
different tapping strengths. Our main result is that the relative contribution of the three different contact types possible between 
tetrahedra - face-to-face (F2F), edge-to-face (E2F), and point contacts - is one variable which stores the preparation history. We 
show this by preparing pairs of packings that differ either in their bulk volume fraction (j) g i 0 bai or in their number of mechanical 
constraints per particle C , where C is determined from all three contact types which each fix a different amount of constraints. For 
the pairs with the same (j)global the number of E2F and point contacts varies with preparation, while the number of F2F contacts 
stays constant. For the iso-constraint packings the relative contribution of all three contact types stays the same. We then perform 
a local analysis of the contact distribution by grouping the particles together according to their individual volume fraction (j)i OC ai 
computed from a Voronoi tessellation. We find that the probability distribution of (j)i OC ai depends only on (j) g i 0 b a u not on C. The 
number of F2F and E2F contacts increases in all experiments with (j)i OC aU the number of point contacts, while always being the 
largest, decreases with (\>i OC ai • However, only the number of F2F contacts can be described by an universal function of (j)i OC ai • 
This behavior differs from spheres and ellipsoids and posses a significant constraint for any mean-field approach to tetrahedra 


packings. 

1 History dependence in granular mat¬ 
ter 

The mechanical properties of a granular sample are known to 
depend on the way the sample was prepared. Examples in¬ 
clude the volume response to sheai^ 1 , the increase of pressure 
with depth in a granular colurm@U, or the density distribution 
insideEl respectively the pressure distribution 7 below a sand- 
pile. Another type of history dependent behavior which is 
sometimes referred to as memory effect can be seen if a sam¬ 
ple is compactified, e.g. by shearing or tapping with a driving 
strength s^, to a specific volume fraction (f) g iobab If then the 
driving strength is changed to so it is found that the subsequent 
evolution of <p g i 0 bai differs for different values of A sim¬ 
ilar effect can be found in granular gases 10 . 

In some sense the term history dependence simply expresses 
the fact that alone does not provide a complete description of 
state of the packing. Given that the mechanical properties of 
the sample originate from forces transmitted between particles 
at contacts it seems natural to extend the description by includ¬ 
ing information on a) the number and spatial structure of the 
contacts and b) the distribution of contact forces. The second 
option is directly connected to the fact that frictional packings 
at finite pressures are hyperstatic 11 14 i.e. there exist many pos¬ 
sible force configurations which fulfill a given boundary con- 
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ditions 15 . The distribution of contact forces has presently only 
be measured in two-dimensional systems of compressive 16 or 
photoelastic 17 discs, the latter experiment finds a clear depen¬ 
dence on the preparation protocol. In principle the preparation 
history could also be stored in distribution of the Voronoi vol¬ 
umes surrounding each particle, but at least for spheres neither 
the volume distribution 18 nor the shape distribution 19 showed 
any signature of the preparation conditions. In this work we 
will study tetrahedra packings both with respect to the number 
and type of contacts and the properties of their Voronoi vol¬ 
umes. 

Several theoretical approaches to sphere packings are based 
on the average number of contacts Z of each particle. For 
example the jamming paradigm! 2021 assumes that mechanical 
properties depend only on the difference A Z between the ac¬ 
tual contact number and the minimum number required to con¬ 
strain all particles (the so called isostatic condition). AZ is 
then assumed to depend on (\) g i 0 bai only. However, for frictional 
particles the formation of contacts is controlled by parameters 
defined on the grain level 22 . This is for example taken into 
account by a statistical mechanics approach to static granular 
media 23 where Z is computed from a mean-field approxima¬ 
tion of (\>iocai • Similarly, the granocentric modeP®D predicts 
the probability distribution of Z in jammed emulsions and disk 
packings using grain level arguments. 

But even the notion that (j) g i 0 bai and Z amount to a com¬ 
plete description of sphere packings has been challenged. The 
value of Z at the Random Close Packing (RCP) density de¬ 
pends e.g. on the numerical protocol used to prepare the sphere 
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3 EXPERIMENTAL SETUP 


packing^. And for frictional disks it has been shown that in 
the cj>global range between RCP and dilatancy onsetf^ mechan¬ 
ical stability requires a fabric tensor anisotropy larger than a 
certain threshold 29 . 


2 Packings of tetrahedra 

History dependence is not limited to sphere packings, recently 
we demonstrated that packings of tetrahedra with the same 
4>global can differ in their values of C and Z depending on how 
the sample was prepared 30 . Such tetrahedra packings have re¬ 
cently been an active field of research with experimental studies 
of their mechanical properties under uniaxial compression 3132 
and their limits of mechanical stability 33 . The jamming of fric¬ 
tionless tetrahedra has been addressed by both simulations 34 41 
and a mean field approach 42 43 . 

Contrary to spheres, tetrahedra have 4 different types of 
contacts: face-to-face (F2F) contacts, which are mechanically 
equivalent to 3 individual point contacts, edge-to-face (E2F) 
contacts (equivalent to two point contacts) and the vertex-to- 
face (V2F) and edge-to-edge (E2E) contacts which we group 
her as point contact. This distinction is important because the 
different contact types fix each a different number of degrees 
of freedom. The total number of constraints C equals therefore 
CiZi where i goes over the three contact types (F2F, E2F, 
point), Zi is the respective contact number and Ci are the num¬ 
ber of constraints per particle fixed at this type of contact. The 
different values of C % can be understood for frictional contacts 
as follows: All contacts impose 3 translational constraints, E2F 
contacts add 2 rotational constraints, F2F contacts prohibit 3 
different rotations. As these constraints are shared between two 
tetrahedra, we obtain the constraint multipliers Cf 2 F = 3.0, 
Ce 2 F = 2.5, and Cy 2 F = Ce 2 E = 1-5. 

Each tetrahedron has 6 degrees of freedom (3 translations 
and 3 rotations), therefore a mechanical stable packing requires 
at least (7 = 6, the so called isostatic case. However, typi¬ 
cal experimental packings of frictional tetrahedra have C val¬ 
ues in the range 12 to 18 and are therefore strongly hyper- 
static^. Contrary claims of isostaticity— are based on the 
use of frictionless constraint multipliers when analyzing exper¬ 
imental i.e. frictional packings. Because of the importance of 
the constraint number for the mechanical properties, we will 
discuss in the remainder of this paper C rather than Z. 


3 Experimental setup 

Experiments were performed with the same polypropylene 
tetrahedra as in reference^ 30 . The particles, figure |l] a) shows 
an example, were made by mould casting, their edge length is 
7 mm and their coefficient of static friction equals 0.8® . 
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Figure 1 Tetrahedra packing, a) Polypropylene particle made by in¬ 
jection moulding, the edge length is 7mm. b) Voronoi Cell around a 
particle, c) Rendering of an interior section of a tetrahedra packing 
after all particles have been identified from an X-ray tomogram. Par¬ 
ticles are colored according to the number of constraints fixed by their 
neighbors. The global volume fraction is 0.54, the average constraint 
number is 14.1. 


Packings were prepared by tapping, starting from a loose ini¬ 
tial configuration which was prepared by first filling the parti¬ 
cles into a smaller inner cylinder (without bottom) and then re¬ 
leasing them into a larger cylindrical container of diameter 10.4 
cm by slowly moving the inner cylinder upwards. Tapping is 
done with an electromagnetic shaker (LDS model V555) which 
is driven by a series of sinusoidal pulses with a frequency / = 
3 Hz. The separation between individual pulses is 0.5 s, the 
amplitude of the sinusoidal motion is chosen such that the peak 
acceleration T corresponds to either 2g, 5 g, or 7g where g is 
the acceleration due to gravity. 

To monitor the evolution of (\) g i 0 bai during tapping, height 
profiles of the packing are measured in regular intervals with 
a laser distance sensor (MicroEpsilon ILD1402) mounted on a 
horizontal translation stage. These readings are then corrected 
for boundary effects by calibrating them with the bulk (pgiobai 
values measured from the tomographic reconstructions. After 
the packings have been prepared, three-dimensional images of 
the geometrical arrangement of particles are obtained by X-ray 
tomography (Nanotom, GE) with a resolution of 100 fim per 
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4 COMPARING GLOBAL VALUES OF ISO-CONSTRAINT AND ISO-VOLUME FRACTION PACKINGS 



Number of taps + 1 


17 


o 


16 - 


O 

0 

o 

Q. 

i_ 

CD 

CL 

(A 


■c3 13 h 

i— 

GO 
C 


11 


- 1 - 1 - 1 - 

) 

2 g 20 ^ 

2 g 1600 

2 g 400 ♦ 

• A 

T 

’ loose ° ^ ^ 0^2 

. O -c. .»- | 

group 1 r 

_i_i_i_ 

^ m 

rUOk ' 

6g 

1 

7 1 Ok 

_i_i_ 


0.46 0.48 0.5 0.52 0.54 0.56 0.58 

Volume fraction 


Figure 2 Compaction of tetrahedra packings for three different tapping 
strengths T = 2g,5g and 7 g. 


voxel. 


Figure 3 History dependent preparation. By choosing the number and 
strength of taps during preparation, we can prepare iso-constraint pairs 
of experiments with different fglobal at similar C (indicated by hor¬ 
izontal arrows) and iso-volume fraction pairs which differ only in C 
not in $ global (vertical arrows). 


3.1 Image analysis 

Particles are identified by a two step algorithm involving a 
cross-correlation with an inscribed sphere and a steepest as¬ 
cend gradient search; the particle detection rates are better than 
99.8 %. After the positions and orientations of all particles are 
known, the number and type of contacts is computed for each 
particle. Figure [T] c) displays a rendering of a cross section 
through the bulk of a sample where the individual tetrahedra 
are color coded according to their number of constraints. Fur¬ 
ther details on these processing steps can be found in the sup¬ 
plementary information of 30 . Similar to reference 30 all particle 
positions are available for download from Dryad®. 

3.2 Preparing packings 

Figure [2] demonstrates the evolution (j) g i 0 bai with the number 
of taps. For a tapping strength T = 2g there is a monotonic 
increase in f g i 0 bai • In contrast, at V = 5g and 7g the volume 
fraction first grows and then decreases and finally grows again. 
In none of the experiments a steady state is reached; in refer- 
ence^t is shown that at least 10 5 taps are needed at T = 2g for 
reaching a plateau in </> g i oha i. 

As already shown in referencd^the number of constraints in 
a tetrahedra packing does not only depend on (j) g i 0 bai but a l so 
on the tapping protocol used to prepare the sample. This allows 
us to prepare two pairs of iso-volume fraction samples i.e. they 
have same the value of f g i 0 bai but which differ in C. Figure [3] 
shows these two pairs, one of them is labeled 7g 10k and 2g 4 , 
the other is 5 g 10k and 2g 1600 . The main number number in¬ 
dicates the T value used during preparation, the superscripts 


denotes the number of taps. 

Our second comparison is between two pairs of packings that 
have approximately the same value of C, but different volume 
fraction <p g i 0 bai • We will refer to them as iso-constraint pairs. 
They are also displayed in figure [5] using the labels 7 g 10k and 
loose (i.e. packing has not been tapped at all); the second pair 
is 5 g 10k and 2 g 20 . 

4 Comparing global values of iso¬ 
constraint and iso-volume fraction 
packings 

The tables in the top right panels of figures [4] and [5] list the 
f global an d C values of all six experiments discussed in this 
paper. The tables also show the contribution of the three differ¬ 
ent contact types towards the total number of constraints. In all 
experiments the overwhelming majority of constraints is con¬ 
tributed by the point and E2F contacts. 

For an evaluation of the relative changes within the pairs an 
estimate for the experimental error is required. As shown in 
the supplements of reference 30 the detection of the total num¬ 
ber of contacts has an approximate statistical error of ± 0.2. 
Unfortunately, we can not derive an error estimate for the algo¬ 
rithm assigning the different contact types. We will therefore 
consider ± 0.2 Ci as a best case approximation of our actual 
errors; this corresponds to db 0.3 for C vo int , ± 0.5 for Ce 2 F , 
and ± 0.6 for Cp 2 F- 

Within these error estimates the relative contribution of the 
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7 MULTIPLICITY OF F2F CONTACTS 


three contact types seems to stay constant for the iso-constraint 
pairs, at least in a first approximation. (Comparing 5 g 10k and 
2 g 20 there might be a small conversion of point to F2F con¬ 
tacts.) For the iso-volume fraction pairs the increase in C is 
only due to an increase in the number of point and E2F con¬ 
tacts, not the F2F. 

The main conclusion of this global analysis is that theoreti¬ 
cal approaches, such as the jamming paradigm, which assume 
that C is determined by (j)global alone will not be sufficient to 
describe packings of frictional tetrahedra. 

We will therefore turn to the local-contact paradigm which 
assumes that the number of contacts a given particle forms will 
only depend on locally defined variables. These can be either 
material parameters, such as the fiction coefficient or the as¬ 
pect ratio a of an ellipsoid. Or parameters which describe the 
neighborhood of the particle, such as volume or shape of its 
Voronoi cell. Here the former is often expressed as the local 
volume fraction fi OC ai which is the particle volume divided by 
the volume of the Voronoi cell. It has recently been shown ex¬ 
perimentally that in packings of spheres^- 7 and ellipsoids^ 22 ! the 
most likely contact number of a particles does only depend on 
f local (and a in the case of ellipsoids), not on fgiobai • Exam¬ 
ples of theories compatible with the local-contact paradigm are 
e.g. the statistical mechanics approach to static granular me¬ 
dia 23 42 or the the granocentric model^ 25 ! 

5 Distribution of local volume fractions 

For both spheres and ellipsoids the probability distribu¬ 
tion of the Voronoi volumes can be rescaled by ( fi OC ai — 
(pgiobai) / v(4>giobai) with cf^giobai) being the standard devia- 
tion of the distribution. The resulting master curve are inde¬ 
pendent of M 

To test the similarity between sphere and tetrahedra packings 
we have computed the Voronoi volume distribution all six ex¬ 
periments; it is displayed in the lower row of panels in figures 

0]and[5l 

6 Contact probabilities depend on local 
volume fraction 

The middle row of panels in figures [4] and [5] displays how the 
contribution of the three different contact types changes with 
4>local • As expected, for all experiments and values of 4>i OC ai 
the number of constraints fixed by point contacts is the largest 
contribution. Interestingly however, the slope dC po i n t / d<pi OC ai 
is negative while it is positive for the E2F and F2F contacts, 
as it is for contacts between spheres or ellipsoids. This agrees 
with the intuitive notion that the closer two tetrahedra get to 
each other, the more likely it is that their flat faces or straight 



Figure 6 Constraints due to F2F contacts in all six experiments. In¬ 
cluded are all bins in fi OC ai which amount to at least 1.7 % of the total 
number of particles. 

edges aligned with each other and the corresponding contact 
changes from a point contact to a higher constrained one. 

At the local level the signature of preparation dependence 
is again clearly visible for the E2F and point contacts; in the 
iso-constraint case their C^fiocai) curves shift vertically while 
preserving their average slope. In the iso-volume fraction case 
one would expect a horizontal shift of the curves. This is how¬ 
ever due to the statistical noise and the rather small slopes only 
evident in the E2F case. 

The only contact type that might behave similar to sphere 47 
and ellipsoid 22 in that it is described by a singular function 
Ci(4>iocai) independent of preparation are the F2F contacts. 
Figure [6] unites the six individual curves. With the possible 
exception of the untapped, loose sample the other curves seem 
indeed to be falling on one master curve. 

7 Multiplicity of F2F contacts 

The different nature of the F2F contacts is not only visible in 
their history independent Cf2F(fiocal ). Figure [7] shows that 
the number of F2F contacts on any given particle are well fit by 
a binomial distribution (with % 2 ^ 10 -3 ); this is not the case 
for the E2F and point contacts. 

This result implies that while the probability of forming an 
F2F contact grows with the proximity between particles, the 
formation itself is an independent random process, lacking any 
geometrical correlation between several neighbors. Such a be¬ 
havior could be rationalized as follows: once an F2F contact is 
formed all three rotational degrees of freedom are completely 
blocked for each particle. Any additional F2F contact depends 
therefore solely on the capability of the incoming particle to 
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align itself such that there is a 180 degree angle between the 
face normal vectors. This is however the same probability that 
lead already to the formation of the first F2F contact. In con¬ 
trast, particles which have already established an E2F or point 
contact retain one or even three rotational degrees of freedom 
to facilitate another contact. 

8 Conclusions 

Because tetrahedra have four flat boundaries, their packings 
do not only posses point contacts similar to sphere packings 
but also edge-to-face and face-to-face contacts which fix larger 
numbers of geometrical constraints than point contacts. The 
number of point and edge-to-face contacts a given tetrahedron 
forms can neither be predicted from the global nor from the lo¬ 
cal volume fraction alone; it does depend also on the prepara¬ 
tion history. In contrast, the formation of face-to-face contacts 
seems to be a stochastic process; its probability is small and 
increases with increasing proximity between particles. While 
our results support the feasibility of a mean-field approach to 
tetrahedra packings, they imply that the jamming paradigm in 
its present form will not be an adequate description. 

We acknowledge computational help by Fabian Schaller. 
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Figure 4 Iso-volume fraction pairs. In all panels the open symbols correspond to the experiment with the higher C. The local analysis in 
the middle row of panels shows how the contribution of the different constraint types depends on the local volume fraction. Nearly the total 
increase in C for the 2 g 400 and 2g 1600 packings is contributed by the E2F and point contacts. At the same time the slope dC/dfiiocai stays 
constant. Ci values are computed for all bins with more than 1.7 % of the totally analyzed particles. The bottom row of panels show the 
probability distribution of 4>i OC ai which suggest little influence of preparation. 
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Figure 5 Iso-constraint pairs. In all panels the open symbols correspond to the experiment with the lower (j) g i 0 bai • The local analysis in the 
middle row of panels shows how the contribution of the different constraint types depends on the local volume fraction. While Ce 2 f versus 
4>iocai display a horizontal shift with increasing (f) g iobai , systematic changes for the point and F2F contacts are more difficult to deduce. Ci 
values are computed for all bins with more than 1.7 % of the total number of analyzed particles. The bottom row of panels show the probability 
distribution of 4>i OC ai which broadens with increasing <fi g iobai- 
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Figure 7 Probability of particles having multiple contacts of the same type, (a) Frequency distribution of the number of face-to-face contacts. 
Straight lines are fits with a binomial distribution resulting in p values of 0.044 (loose), 0.054 (2 g 400 ), and 0.064 (7 g 10k ). In contrast, the 
multiplicity of the E2F (panel b) and point contacts (panel c) is not described by a binomial distribution. 
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